Enhance synchronizability via age-based coupling 
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In this brief report, we study the synchronization of growing scale- free networks. An asymmetrical 
age-based coupling method is proposed with only one free parameter a. Although the coupling 
matrix is asymmetric, our coupling method could guarantee that all the eigenvalues are non-negative 
reals. The eigneratio R will approach to 1 in the large limit of a. 
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One of the main goals in studying of network science 
is to understand the relation between network struc- 
ture and dynamical processes performed upon [H, 0] . A 
typical collective dynamic on networked system is syn- 
chronization, where all the participants behave alike, 
even exactly the same. This phenomenon exists every- 
where from physics to biology 0], and has been ob- 
served for hundreds of years. With the partial knowl- 
edge of relations between network structure and its syn- 
chronizability 0, H, H, 0, HI, scientists have proposed 
many methods to enhance the network synchronizability 
H El HA El, El, 0, El El El, EH . Generally speaking, 
these methods can be divided into two classes, one is to 
modify the network structure I9l lid-Jill , the other is to 
regulate the coupling pattern [l3.ll3l1l4l flBI . El El, El] • 
In the former class, networks are modified either to 
shorten the average distance [l(| or to eliminate the max- 
imal betweenness d E3 . In the later case, the network 
structure is kept unchanged, while the coupling matrix is 
elaborately designed (often asymmetric) to improve the 
synchronizability 



led ott en asymmetric) to 

El ES Q El, El, [13, El- 



The first coupling pattern other than the symmetric 
case was proposed by Motter-Zhou-Kurths [H, El, E3 
(MZK coupling pattern), in which the coupling strength 
a node i receives from its neighbors is inverse to fef with 
ki the degree of i. The coupling pattern can sharply en- 
hance the network synchronizability, with /3 = 1 the op- 
timal case. After this pioneer work, many coupling pat- 
terns [IH, El El, [H| have been presented to further en- 
hance the network synchronizability. In Ref . [151 ] , Hwang 
et al. presented a coupling method taking into account 
the age of nodes, which makes the network even more 
synchronizable than the optimal case of MZK coupling 
pattern. In this pattern, each node receives coupling 
signals from its neighbors, with each receiving coupling 
strength taking one of the two values: if the neighbor is 
older, the coupling strength takes the larger value, oth- 
erwise it takes the smaller one. To separate the differ- 
ent coupling situations (i.e. from older to younger and 
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from younger to older) by using two different coupling 
strengths is the simplest way one can image. However, 
since each node has its own age, a coupling method tak- 
ing into account the age difference between each pair of 
coupled nodes may further enhance the synchronizability. 
Moreover, the coupling matrix in Ref. [151 ] has complex 
eigenvalues, leading to a complicated analysis. An elab- 
orately designed method, as shown in this brief report, 
could guarantee that all the eigenvalues are nonnegative 
reals, thus one can easily predict the synchronizability 
of underlying network by considering the real eigneratio 
only. 

In a dynamical network, each node represents an os- 
cillator and the edges represent the couplings between 
nodes. For a network of N linearly coupled identical os- 
cillators, the dynamical equation of each oscillator can 
be written as 



N 

F(x*)-a^G y -H(x-'), 

3=1 



1,2, ...,N, 



(1) 



where x l = F(x') governs the essential dynamics of the 
ith oscillator, H(x J ) the output function, a the coupling 
strength, and Gtj an element of the N x N coupling ma- 
trix G. To guarantee the synchronization manifold an in- 
variant manifold, the matrix G should has zero row-sum. 
The collective dynamic starts from a disorder initial con- 
figuration, under suitable conditions, the couplings will 
make the oscillators' states nearer and nearer. Eventu- 
ally, all the individuals oscillate together, and synchro- 
nization phenomenon emerges. 

In the simplest symmetric way, the coupling matrix G 
has the same form as the Laplacian matrix L, that is, 



Gtj = Lij where 



ki for i — j 
-1 for j S Aj 
otherwise 



(2) 



Here Aj is the set of z's neighbors. Because of the sym- 
metry and the positive semidefinite of L, all its eigen- 
values are nonnegative reals and the smallest eigenvalue 
Ao is always zero, for the rows of L have zero sum. 
And if the network is connected, there is only one zero 
eigenvalue. Thus, the eigenvalues can be ranked as 
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Ao < Ai < A2 < ... < Aat_i. When the stability zone 
is bounded, according to the criteria of master stability 
function [HJ Hfj (see also the unbounded case [2ll [22l|). 
the network synchronizability can be measured by the 
eigenratio R = Ajv-i/Ai: The smaller it is the better the 
network synchronizability and vice versa. 

The couplings between nodes are not limited to the 
symmetric mode, however, generally, the eigenratio of an 
asymmetric coupling matrix is complex (e.g. the eigen- 
ratio in Ref. Therefore, in order to ensure the 
network having strong synchronizability, not only the ra- 
tio of the real part should be taken into account, but 
also the imaginary part must be guaranteed as small as 
possible. In Ref. [la], the simulation result indicated 
that although the ratio of the real part is the smallest, 
at the same time the imaginary part is the largest. To 
overcome this blemish and give further enhancement of 
synchronizability, we bring forward a coupling pattern 
in which the coupling strength between two connected 
nodes is the function of their age difference. The age of 
node i is signed by the time it enters into the network, 
thus smaller i corresponds to older age. The coupling 
matrix proposed here is: 



1 



Gi 



Si 





for i = j 

for j G Ai 
otherwise 



(3) 



where Si — X)jga e~^^~ 1 ' is the normalization factor. 
In this coupling pattern, the case of a = degenerates to 
the optimal case of MZK coupling pattern. When a > 0, 
the old nodes have stronger influence than the younger 
ones; while for a < 0, younger nodes are more influential. 

It can be proved that although the coupling between 
nodes is asymmetric, all the eigenvalues of matrix G are 
reals. Note that, the coupling matrix defined in (3) can 
be written as 



G = DL' 



(4) 



where 



D = diag(e*7Si, e 4Q /5 2 , e» a /S 3 , e*» a /S N ) (5) 

is a diagonal matrix, and L' = [L'^) is a symmetric zero 
row-sum matrix, whose off-diagonal elements are 



L', 



(G) 



From the identity 

det(DL' - XI) =det(DH'Di - XI) (7) 

valid for any A, we have that the spectrum of eigenvalues 
of matrix G is equal to the spectrum of a symmetric 
matrix defined as 



H = D^L'D'- 



(8) 



As a result, although the coupling matrix G is asymme- 
try, the eigenvalues of matrix G are all nonnegative reals 




FIG. 1: (Color online) The eigenratio R vs. a in BA networks 
with average degree (k) — 6. The inset displays the details 
for the interval a £ [—4,1]. Each data point is obtained 
by averaging over 50 different network configurations. The 
eigenratio R goes to 1 in the large limit of a. 



and the smallest eigenvalue is always zero. Therefore, 
different from the complicated case in Ref. [l5| , the syn- 
chronizability based on the present coupling pattern can 
be measured directly by the real eigenratio R. 

In Fig. 1, we report the changes of eigenratio R with 
the parameter a in BA networks [23| at different sizes. 
One can easily conclude from Fig. 1 that with the in- 
creasing of a the eigenration decreases sharply, no mat- 
ter what the network size is. It is shown that in growing 
networks, if the couplings from older nodes are stronger 
than the reverse, the network will get better synchroniz- 
ability. Otherwise, if the coupling from younger to older 
ones is strengthened (see the cases of a < in the inset), 
the system becomes very hard to synchronize. When a 
goes to infinite, the eigenratio will converge to 1, which is 
the possibly smallest eigenratio corresponding to the best 
synchronizability [241 ] . Actually, in the case a — > +00, 
each node is coupled by its oldest neighbor, while the 
oldest node in the network is uncoupled. Thus, the cou- 
pling matrix (whose rows are sorted by the descending 
order of ages) becomes a lower triangular matrix with all 
the diagonal elements are 1 except the first one, Gu, be- 
ing equal to zero. Therefore, all the non-zero eigenvalues 
are one. 

Although there exists various methods to design a 
coupling pat tern having optimal synchronizability (i.e. 
R = 1) [241 1 f° r growing networks, using the age of 
each node is a simple and feasible way since to know 
any other measures of nodes may cost much for huge- 
size system, and this age-based coupling can guarantee 
the connectivity of the whole network. Mathematically 
speaking, the synchronizability here is a measure on the 
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stability of invariant synchronization manifold. We call a 
synchronization manifold is stable if the dynamical sys- 
tem can automatically return to this manifold after a 
perturbation. A network G has better synchronizability 
than another network G' means any collective dynamics 
with identical oscillators upon G' having a stable syn- 
chronization manifold must have a stable synchroniza- 
tion manifold for G, while there exists certain dynamics 
having stable synchronization manifold for G, but not for 
G'. However, better synchronizability does not guarantee 
a shorter converging time from disorder initial configu- 
ration to synchronized state. Actually, Nishikawa and 
Motter [2J] found that the synchronizing process may 
take longer time in the optimal network with R = 1 (see 



also a similar conclusion for non- identical oscillators [H|). 
Based on the current coupling pattern, one can obtain his 
acceptable trade-off between synchronizability and con- 
verging time by tuning the parameter a. Moreover, com- 
paring with the pioneer work by Hwang et al. [la ], our 
coupling method can achieve even smaller R, and does 
not need to deal with the complicated and boring anal- 
ysis on complex eigenratio. Instead, our elaborately de- 
signed coupling pattern can guarantee the eigenratio a 
real number. 
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